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Abstract 

Making a connection between observations of cosmological correlation functions and those cal- 
culated from theories of the early universe requires that these quantities are conserved through the 
periods of the universe which we do not understand. In this paper, the results of [0810.2831] are 
extended to show that tree-approximation correlation functions of Heisenberg picture operators for 
the reduced spatial metric are constant outside the horizon during local thermal equilibrium with 
no non-zero conserved quantum numbers. 
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I. INTRODUCTION 



If observations of the cosmic microwave background or large scale structure are to tell us 
anything about the early universe, we need to understand how the correlation functions that 
we calculate from theories of inflation are related to the correlation functions that are present 
during the radiation-dominated era at a temperature low enough such that we understand 
the contents of the universe. This presents a problem because there are several eras in the 
early universe which we do not understand including dark matter decoupling and baryon 
and lepton synthesis. 

In [l] it was shown that under essentially all conditions the non-linear classical equa- 
tions governing matter and gravitation in cosmology have adiabatic solutions in which, far 
outside the horizon, in a suitable gauge, the reduced spatial metric <7y(x, i)/a 2 {t) becomes 
a time-independent function Qij(jc), and all perturbations to the other metric components 
and to all matter variables vanish. Furthermore, it was shown that the adiabatic solution 
is attractive in the case of local thermal equilibrium with no non-zero conserved quantum 
numbers. However, this is not enough to conclude the the correlation functions for ^(x, t) 
become time-independent at late time. In calculating non-gaussian correlations of the re- 
duced metric, one generically encounters loop graphs which receive contributions from fields 
with arbitrarily large wave numbers, invalidating the expansion in powers of a" 1 as applied 
to the Heisenberg picture interacting fields [2]. In fact, this problem arises even for classical 
fields, where the initial conditions for the non-linear field equations are taken to be stochas- 
tic variables. Furthermore, the prescription for calculating correlation functions even at the 
leading order in perturbation theory introduces an implicit time-dependence which must be 
shown to vanish in the appropriate limit if correlation functions truly become constant. Here 
we apply the results of [3( to show that in the tree approximation the correlation functions 
of Heisenberg picture operators for the reduced spatial metric are indeed constant outside 
the horizon during local thermal equilibrium with no non-zero conserved quantum numbers. 

II. EXPLICIT ADIABATIC SOLUTION 

Here we review the results of [l| which found a solution of the non-linear classical grav- 
itational field equations for which the reduced metric approaches a constant at late times. 
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Most of this section simply recalls the arguments given in however we will slightly extend 
these results by calculating the subleading terms in the momentum potential Z7(x, t) which 
we will need in section IIHI 

n 

We will use the ADM parametrization |4J for the metric: 

g 00 = -N 2 + gyrfW , g 0l = 9ij N j = N, , 

g 00 = —TV -2 , g° l = N' l /N 2 , g ij = (3) g ij - N' l N j /N 2 , (1) 

where is the reciprocal of the 3 x 3- matrix g^. It will be convenient also to write 

gij(x,t) = a 2 (t)g i:j (x,t) , (2) 

where a(t) is the Robertson- Walker scale factor appearing in the unperturbed solution. If 
we take units with 8ttG = 1 the Lagrangian takes the form 

L = ^Jd 3 x ^-Detg {-R^} + L m 

= °- J d'xN^g [-argufy + C)C\ - {C\f + 2C m ] , (3) 

and the gravitational field equations are 

V, {C) - ^.C\) = -NTf , (4) 

- a- 2 g ij Rij - C)C\ + (C\) 2 = 2N 2 T 00 , (5) 

Rij - C k k C ij + 20*0^ + N- 1 (- Cij + C\W 3 N k + C k 3 WiN k 

+ N k V k C l3 + ViVjN) = -Ty + l -a 2 ~ gi] T\ , (6) 

where <? u (x, t) is the reciprocal of the matrix <7y (x, t); -Rjj(x, t) is the three-dimensional Ricci 
tensor for the metric g^{x,t); and CUx,t) is the extrinsic curvature of the surfaces of fixed 
time 



1 a = a 2 g lk C k j , C ki = 



2aag ki + a 2 g ki - VfciV; - ViN k 



(7) 
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where V« is the three-dimensional covariant derivative calculated with the three- metric g^. 



Following jl], we will look for a solution in which g^ and goo + 1 are small perturbations 
for large a(t), of order l/a 2 (t). We will define space coordinates for which N l = 0, so g^ = 0. 
We then have 

(■'., = H6) + Cj (8) 



where £*■ is, like ^ and 5N, a quantity whose leading term is of order a 2 : 

ft = ^ fe - 2ff<WV0y] + 0(a- 4 ) . (9) 
The total energy momentum tensor for this solution takes the form 

Tij = a 2 pg i:j + 8T i:j , T i0 = 5T i0 , T 00 = p + 5T 00 , (10) 

where the background energy density and pressure are given by p = 3H 2 , p = —2H — 3H 2 , 
and the quantities STij/a 2 , ST i0 , and 5T 00 are all of order a~ 2 . The gravitational field 
equations now read 

^ (ft - = ST 0j + 0(a- 4 ) , (11) 

- 12H 2 5N + 25T 00 = -a'^R^ + 4H£\ + 0(a" 4 ) , (12) 



ft- + 3#ft + #<^£\ + (3# 2 - if^MJV = a~V fc i?, 
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+ a,- 2 g ik 6T kj - -a- 2 5y g kl 5T kj + -5)5T 00 + 0(a~ 4 ) . (13) 

To order or 2 adiabatic solutions have no vorticity, so we can introduce a momentum 
potential U, of order a -2 , such that 

5T i0 = d l U. (14) 
The equation of momentum conservation then reads 

a~ 2 g ik Vi5T kj = 2Hd j 5N + dj(U + 3HU) + 0(a~ 4 ) . (15) 
We can always write 

6T ij = a 2 (g ij 6p + n ij ), (16) 

where is a 3-tensor of order a~ 2 representing anisotropic inertia, with g^Hij = 0. For 
models with vanishing anisotropic inertia, as in the case of local thermal and chemical 
equilibrium, we have 



g lk Rkj — -^ % j9 kl Rki 



+ 0(a" 4 ), (17) 



where 



S , ; = r, + 2^- ( 18 ) 
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The general solution of ([IT]) is 



5*(x,t) 



£ ifc (x)^(x) - -^"(x^x) 



a»(t) 



a(t') rft' 



+ 



^( x ) 
a»(t) 



0(a" 



(19) 



where T is any fixed time, ftj(x) is the value of ^(x, t) at that time, 7£jj(x) is the Ricci ten- 
sor calculated from the 3-metric £?y(x); and -B l -(x) is some function of x (and T), appearing 
in the solution of the homogeneous equation corresponding to (1171) . 

Next, we need to solve for the metric. From (JUJ), (Il2p. and (ITS]) we have 



tl 



2a 2 H 



T -g ij g kl R k i + g ij X + 0(a 



where X = 0(a 2 ) depends on the matter perturbations 

X = --r- (sTqq - 2(3H 2 + H)5N^J -U[l 



(20) 



(21) 



Under a shift t — > t + e(x, t) in the time coordinate, with e of order a~ 2 (and a corresponding 
transformation x % — > x % + £/ u J dt a~ 2 de/dx^ to keep iVj = 0), the quantity X undergoes the 
transformation 

X -+ X + 2^ (eif ) + . . . , (22) 

so we can choose e to make X = 0. This choice does not fix the gauge completely, however. 
There is a residual gauge freedom which preserves both and X = which takes the form 



t -H + r(x)/H(t) 



x i -> + £ ij (x 



<9r(x) 



(23) 



^ J a 2 (t)H(t) ' 
where r is an arbitrary function only of x. 

With the gauge choice X = we can solve for g^ without solving for the matter variables 
to which the metric is coupled. This solution takes the form 



<7i 3 -(x,i) =Gij(x) + 2 



72y(x) - -<fc(x)0«(x)72«(x) 



o^t 7 ) 



a(t") (it'' 



+ ^ l ,(x)^(x)^(x)^ < -^ +/ 



iP(t') dt' 



(t')a 3 (f) 



a(t") dt'' 



xW H x )ft w x / — ^4 — 
2 ^ 1 J " l } J T a 2 (t')H(t' 



+ 2<^,(x) J B fe J (x) y +20 y (x)fl*,(x) 
+ 0(a- 4 (t)) , 



f H 2 (t')dt' 
T a?(t')H(t') 



(24) 
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where T is again any fixed time, and ^(x) and IZij (x) are the values of and the associated 
Ricci tensor at that time. 

Next we need to solve for the remaining metric component g>oo = —N 2 . From the perfect- 
fluid form of the energy momentum tensor, we can write 

5T 00 = 5p + QH 2 5N + 0(a- 4 ) (25) 

so that ([12} and give 

25 p = -a^gVRij + 4HZ k k - QHU + 0(a~ 4 ) , (26) 
while fl25|) and f l2T|) with the gauge condition X = give 



= — r- (£p - 2H5N) - U I 1 + — ^- ) + 0(cT 4 ) . (27) 



With vanishing anisotropic inertia (TT5l) gives 

= -5p + 2H5N + 3HU + U + 0(a- A ) . (28) 

We now have three relations for the four quantities 5N, Sp, Sp, and U , so we need one 
additional relation to calculate all four. 

In local thermal equilibrium with no non-zero conserved quantum numbers the pressure 
and energy density in any gauge are functions only of the temperature, so 

Combining (l26l)-([29l) we arrive at a differential equation for the momentum potential U 



U + U ( ^ - 5- ) = - ( 2 + I {-—J'Rii + 2HE* 

H HI \ 3HH J V 2a 2y ' 

H \ ( 1 H 



-( i+ o4) (-h + §f 

( 4 ^+^H^ + °( a " 4 ), (30) 



3H J a 

where T may be taken as any time during the period of thermal equilibrium, most conve- 
niently at its beginning, and Gij and TZ^ are the reduced metric gij/a 2 and the associated 



Ricci tensor at that time. The solution of this equation is 



-. W , W ||/>(g)( 1 + ^) 

+ 0(a~ 4 ), (31) 



where /(x) is an arbitrary function of position. We can slightly simplify this expression by 
using the relation 



d ( H(t') \ 1 3H 2 (t') H(t')H(t' 



dt' \a 3 (t')H(t')J a 3 (t') a 3 H(t>) a%t')H 2 (t') 
so the fourth line of ( T3TT) becomes 



2H{T) -B^m+^B^ 



(32) 



3a 3 {T)H{T) lV 'H(t) 3a 3 (t)~ 
Vfx^ (t) T dt> fl+3^ ( ^ f33^ 

-^^mu^m \ 1+3 im)- (33) 

The first term can be absorbed into a redefinition of the arbitrary function /(x). The leading 
term in (13TI) . representing a solution of the homogeneous equation corresponding to ( )30|) . is 
of zeroth order in a -1 , and so does not become small for large a. However, this term can be 
removed by the residual gauge transformation ( 1231) under which 



C/(x,t)^C/(x,t)+^^jr(x). (34) 
By choosing r(x) to have the value — /(x)/2, we can cancel the first term in (|3~T1) . After this 
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gauge transformation, U takes the form 



6H(t')H{t>) 



3a3(t) 1W 3 4V ' H{t) J T a?{t>) \ H(t>) 
+ 0(a" 4 ) , (35) 

The leading term of the remainder is of order a -2 and thus vanishes for large a. 

The remaining perturbations 5goo = —25N and 5p are algebraically related to U by ( 126|) 
and (12 7p. which give 

5p(x,t) =-3H(t)U(x,t) 

+ i^(x)^,(x) (~ + || "(*')) + 0(0 (36) 



and 



2H(t)6N(x,t) = ( ^4 I ^(x,t) 



^(0 

+ ^(x)^(x) (--^ + §| /W)) + O(a-) • (37) 

Both are of order a -2 , so for local thermal equilibrium with no non-zero conserved quantum 
numbers, the adiabatic solution of the non-linear classical field equations is attractive for 
large a. 



III. TREE-APPROXIMATION CORRELATION FUNCTIONS 

The results of which we reviewed in the previous section, show that the non-linear 
classical field equations have an adiabatic solution in which the reduced spatial metric is 
constant outside the horizon and all perturbations to the other metric components and 
the matter variables vanish, and this solution is attractive during local thermal equilibrium 
with no non-zero conserved quantum numbers. However, this is not enough to conclude 
that correlation functions of the reduced spatial metric are conserved outside the horizon. 
In calculating correlation functions of Heisenberg picture interacting fields, one generally 
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encounters loop graphs which receive contributions from fluctuations with arbitrarily small 
wave number, invalidating the expansion in powers of a -1 . In this section we will apply 
the general theorem of [3] to show that in the tree approximation, correlation functions of 
the reduced spatial metric are constant outside the horizon during local thermal equilibrium 
with no non-zero conserved quantum numbers. This section will closely follow the analysis 
of 2! which showed that tree approximation correlations are constant during single field 
inflation. 

The generating function for correlation functions of at a time t\ is 



exp{W[J,ti]} = (0,in 



exp 



d*x~g?{xM)J ij {x) 



0,in 



(38) 



where gfj (x, t) is the Heisenberg-picture quantum mechanical operator corresponding to 
(jij(x.,t). Correlation functions for <?y(x, t) are calculated by the formula 

Qn 



0,in l&j (x, ti)<7 w (y,ti 



0,in) 



ex V {W[J,h]} 



(39) 



j=o 



_a/^(x)a/ w (x) 

In order to calculate the generating function W in the tree approximation, we need to 
construct complex c-number metric fields gij(x,t) together with a complex auxiliary field 
N(pc,t), satisfying the following constraints: 

(A) The fields satisfy the Euler-Lagrange equations. For this case, they are (@|-(|6]). 

(B) The fields g^ satisfy constraints at the time t\ when the correlation functions are 
evaluated 



Im^(x,ti) = , 
\ $7<iCMi) 



■J* 3 (x). 



(40) 
(41) 



(C) gij satisfies a positive frequency constraint at time t — > —00, so that it behaves as 
a superposition of terms proportional to exp(— iut), with u various positive frequencies. 
(This final constraint would be modified if one considered a non-trivial initial state for the 
cosmo logical fluctuations.) 

With gij(pc,t) and N(pc,t) calculated subject to these constraints, the contribution of 
connected tree graphs to the generating function is given by 

W[J,h] tiee = f lmL[g(t),l(t),t] dt+ fcPx^Wgifah). (42) 

J —CO J 



These constraints give the functions g^ an implicit dependence on the time t% when the 
correlations are evaluated, so in order to conclude that the correlation functions are truly 
time- independent, we must show that the constraints become independent of t\ for suffi- 
ciently large a(ii). Furthermore, we must show that the time integral appearing in ( HSj) 
converges for sufficiently large a(t\). 

For large a(ti), the constraint (T4"0"|) provides the ti-independent condition that leading 
term Qij(x.) in the solution for the reduced spatial metric fT24"j) must be real for all x. The 
associated Ricci tensor IZij (x) is then also real, so all of the terms of order a" 2 appearing in 
ff24"|) are real. As a result 

Im^(x,t)=2^(x)Im^.(x) / — + 2ft J (x)ImB\(x) / 

Jt a [t ') J T a 3 (t')H(t') 

+ 0(a-\t)), (43) 

and the leading terms in Img(x, ti) are of order a~ 3 . 

We will next check the convergence of the time integral appearing in (T4"2|) for large t\, 
returning below to the functional derivative appearing in the constraint (|4ip . In order to 
do this, we need to say something about the matter Lagrangian L m appearing in ([3]). As 
long as comoving entropy density is conserved, which is guaranteed by the condition of local 
thermal equilibrium with negligible chemical potentials, apart from some total derivative 
terms the value of matter Lagrangian is given simply by 5rtZ] 

L m = J d 3 x ^-Detgp = a 3 J d 3 x [p + 5p) , (44) 

where p is the pressure. 

By examining fl24"|) . fl35|) . f l29|) . fl36l) . and (1371) we can see that each of gij, U, 5p, and 
5N are quantities of order a~ 2 with leading imaginary part of order a~ 3 . Any second order 
function of these quantities will thus have leading imaginary part a~ 2 x a~ 3 . 

We will first consider the terms in the Lagrangian ([3]) containing either or 1 space or 
time derivative and which are of zeroth order in 5N and 5p. These terms are 



= Yj d3x ^ [- 12H2 2H 9 ij lij\ , (45) 
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where we have used the fact p = —2H — 3H 2 . As was pointed out in [8], this is a total time 
derivative 



L [g,g,t] 



-2j[ /' /r/ :; r v A/) . 



(46) 



It was shown in the appendix of [2] that a total time derivative in the Lagrangian has no 
effect on the correlation functions, so we can ignore these terms in what follows. Next, we 
will consider the terms of zeroth order in and which are of first order in SN or Sp. These 
terms are 

,3 



Li [g, M = y / £x \T$ [<°H 2 5N + 2pSN + 2Sp] 



— I d 6 x ^g -4H5N + 2Sp 



(47) 



We can add to the Lagrangian a term 



ALi [g, g, t] 



3 'd 3 x^U 



3HU + U + -U~g ij tij 



d 

Tt 1 " 



(48) 



which is a total derivative so it does not alter the correlation functions. Combining this with 
L\, we have 



(Li + ALi) [g,g,t] = y / d 3 x ^g 



a 

T 



-AH5N + 25p - 6FC7 — 2Z7 — U~g l3 g ij 
d 3 x^ (U~g l ^ giJ + 0{a- A )) , 



(49) 



where we have used (I28p to simplify the second line, and the 0(a~ 4 ) term is a second order 
function of U, 5N, and Sp. The leading imaginary part of L\ + AL\ is thus of order a~ 2 . 
Next, we will consider the term in L containing the spatial curvature. By examining fl2] 
and (I37p . we can see that the leading imaginary part of the term 

(50) 



L s [g,g,t] = -— j d 3 xN^a 2 ~g ij R i3 

is of order a~ 2 . All of the remaining terms in the Lagrangian are given by a 3 times some 
second order function of c/ij, U, Sp, and SN, plus terms of higher order in a~ l . As a result, 
all of the remaining terms have a leading imaginary part of order a~ 2 . We therefore conclude 
that the time integral appearing in (|42p converges to a finite limit for large t 1 as J 4l a~ 2 (t) dt. 
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Finally, we will consider the functional derivative appearing in the constraint (|4T|) . For 
this purpose, we will add to the Lagrangian the total derivatives — L and ALi identified 
above in (1461) and (l4"gj) . so we have 
5{L-L + AL 1 ) [g{t),l{t),t] 

= ^f^ g^x, t) {-2H{t)5\ + e fc (x, t) - ^(x, t)) 
_ 5L [g{t)M)A + SALx [g(t),t(t),t] 

= ^f^j'lx, t) (-U(x, + f' fe (x, f) - <^(x, + 0(a~ 4 )) . 

Recalling the definition ( TT8|) . we can rewrite this as 

g(L-Lp + ALi) [fl(t),5(t),t] 
5gij(x,t) 

V(x, t) fe\(x, t) - <^H< (x, t) + 0(a- 4 )l . (52) 



2 

The terms of order a -2 appearing in ( JT9l) are real, so the leading imaginary part of H*- is 
given by 

ImS^x^) = a-^ijlmB^x) +0 (a~ 4 (t)) . (53) 

As a result, we have 

5(L-L + ALi) [g(t),§(t),t] 



Im 



8gij(x,t) 



= ^f^(x) Im (^ fc (x) - ^S',(x)) + O {a-\t)) , (54) 

and so the constraint ( l4"Tj) becomes independent of ti for large a(ti). 

We have shown that during local thermal equilibrium with no non-zero conserved quan- 
tum numbers the tree-approximation generating function Wtreef^i] converges to a t\- 
independent function for large a(ti), and so the correlation functions calculated in the tree 
approximation are constant far outside the horizon. 



IV. CONCLUSIONS 



in 



We have used the results of [1| and [3] to show that the tree- approximation correlation 
functions of the reduced spatial metric are constant outside the horizon during local 
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thermal equilibrium with no non-zero conserved quantum numbers. This result shows that 
if we are able to calculate correlation functions from a particular theory of inflation and 
follow the evolution up until a phase of local thermal equilibrium, then we may compare the 
predictions of the theory with observations made today, even without an understanding of 
several phases of the early universe which follow inflation, as long as loop corrections can 
be neglected. It has been shown that loop corrections are much smaller than the leading 
contributions in typical theories of inflation jsljH], and so the tree approximation is sufficient 
for making a connection with observations. 

The results presented here are exactly what would be expected from the intuition that 
the evolution in the tree approximation should correspond to the classical solution. On the 
other hand, the explicit proof given here guarantees the constancy of cosmological correlation 
functions without reference to the classicality of the relevant quantities. 
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